1 



Finite size scaling analysis of compact QED 

G. Arnold, Th. Lippert, K. Schilling 

Department of Physics, University of Wuppertal, D-42097 Wuppertal, Germany 
Th. Neuhaus 

Department of Physics, University of Aachen, D-52056 Aachen, Germany 



o 
o 
o 

(N 
> 

o 



> 

00 

o 



o 
o 



X 



We describe results of a high-statistics finite size scaling analysis of 4d compact U(l) lattice gauge theory with 
Wilson action at the phase transition point. Using a multicanonical hybrid Monte Carlo algorithm we generate 
data samples with more than 150 tunneling events between the metastable states of the system, on lattice sizes 
up to 18*. We performed a first analysis within the Borgs-Kotecky finite size scaling scheme. As a result, we 
report evidence for a first-order phase transition with a plaquette energy gap, G = 0.02667(20), at a transition 
coupling, /3t = 1.011128(11). 



1. INTRODUCTION 

The determination of the order of phase tran- 
sitions is of high importance for lattice field theo- 
ries. For it requires higher than first-order phase 
transitions to make contact between lattice and 
continuum physics. 

In compact quantum electrodynamics (QED) 
with Wilson action the order of the transition be- 
tween the confined and Coulomb phases has been 
under debate for many years On the scale of 
accessible lattice sizes the correlation length is 
large, but the system exhibits a definite two-peak 
structure. 

In a recent high statistics investigation [Q the 
latent heat appeared to decrease with the lattice 
size L, with a critical exponent being neither 
0.25 (first order) nor 0.5 (trivially second order), 
with significant subasymptotic contributions at 
the studied values oi L < 12. These findings allow 
for two scenarios: 

1. The observed double peak structure is a 
finite size effect and vanishes in the ther- 
modynamic limit. The signature of a 
second-order phase transition would even- 
tually emerge at some Lq ^ 12. 

2. The phase transition is weakly first-order, 
i.e., the correlation length ^ remains finite, 
yet large in terms of lattice extensions ac- 
cessible today; this would fake, on small lat- 
tices, the signature of a second-order transi- 



tion, since the true value of the latent heat 
would be revealed only in the regime L > ^ 



For a class of spin models with strong first- 
order phase transitions finite size scaling the- 
ory has become amenable to quantitive stud- 
ies through the work of Borgs and Koteckyj^ 
||. There the finite volume partition function 
at temperature /3 in finite volumes with periodic 
boundary conditions (neglecting interfacial con- 
tributions) has the remarkably simple form 



Z = \e-^f^(f>'> + g-V/2(/3)+ln(X)l 



(1) 



The functions /i(/3) and f2{f3) denote bulk free 
energy densities in the two coexisting phases 1 
and 2. X denotes the asymmetry parameter 
which is nothing but the relative phase weight 
in the probability distribution P{E). 

In Ref. H] the Borgs-Kotecky ansatz has been 
extended in a heuristic manner to the 3d 3- 
state Potts model undergoing a weakly first-order 
phase transition; in this instance detailed con- 
sistency checks have been carried out in order 
to verify the viability of the Borgs-Kotecky ap- 
proach. Motivated by this success we shall apply 
in the following this ansatz to compact QED lat- 
tice gauge theory. 



2 



2. SIMULATION DETAILS 

We consider 4d pure U(l) gauge theory with 
standard Wilson action 



S=-f3 



cos(6'^^(n)), 



(2) 



where P denotes the Wilson coupling and 9^i,{n) 
the plaquette angle. We use a cubic lattice of vol- 
ume V = with periodic boundary conditions. 

We have implemented three different algo- 
rithms for generating the U(l) gauge field config- 
urations: (a) a local Metropolis (MRS), updating 
each link separately, (b) a global hybrid Monte 
Carlo algorithm (HMC) and (c) a combination of 
the multicanonical and the hybrid Monte Carlo 
algorithm (MHMC). For details we refer to Refs. 
^,01 . Each update of the complete lattice is fol- 
lowed by 3 reflection steps |^ to reduce correla- 
tion of successive configurations. The number of 
generated configurations at each lattice size L is 
> 2.5 X 10^. We additionally measure the number 
of tunneling events (flips) as control parameter for 
the mobility of the algorithms. 

Our simulation parameters and the statistics 
achieved are listed in Table ^ As the runs differ- 
ing by algorithm, HMC parameters or by weight 
function are independent we do not recombine 
them into one multihistogram. Rather we treat 
them separately. 



3. MEASUREMENTS 

Based on the plaquette operator 

E cos(0^.(n)). 



n,i/>/j, 

we consider the following cumulants: 
a(/3,L) = &p^v({E^)-{Ef), 



U2iP,L) = 1- 



{EY 



UiW,L) = 1(1- i:^ 

3 I {E^)\ 



(3) 

(4) 
(5) 

(6) 



The locations {fic,C^v,raax) and {Pc,Ui^mi,i) of 
their extrema are determined by reweighting the 



Table 1 

Number of configurations generated on different 
lattice sizes. In case of HMC generated data the 
time step At and trajectory length Nmd is speci- 
fied. The MHMC runs with equal parameter set- 
tings differ by the choice of weight parameters. 



1 1 i " i i m. n. i i —r-- , 


6 


1.001600 


HMC 


. 120 


10 


2500000 


3581 


8 


1.007370 


HMC 


.093 


2 


1250000 


276 




1.007370 


HMC 


.093 


4 


1250000 






1.007370 


HMC 


.093 


6 


1250000 


707 




1.007370 


HMC 


.093 


8 


1250000 


856 




1.007370 


HMC 


.093 


9 


1250000 


913 




1.007370 


HMC 


.093 


10 


1250000 


932 




1.007370 


HMC 


.093 


11 


1250000 


935 




1.007370 


HMC 


.093 


12 


1250000 


921 




1.007370 


HMC 


.093 


13 


1250000 


935 




1.007370 


HMC 


.093 


14 


1250000 


859 




1.007370 


HMC 


.093 


16 


1250000 


862 




1.007370 


HMC 


.093 


13 


1440000 


1265 


10 


1.009300 


HMC 


.071 


9 


1000000 


294 




1.009300 


HMC 


.071 


11 


1000000 


350 




1.009300 


HMC 


.071 


15 


1000000 


351 




1.009300 


HMC 


.071 


17 


1000000 


344 




1.009300 


HMC 


.071 


19 


1000000 


340 


12 


1.010143 


MRS 






3617000 


571 






MHMC 


.060 


20 


1303500 


407 


14 


1.010598 


MRS 






3900000 


215 




1.010568 


MHMC 


.050 


24 


825200 


186 


16 


1.010753 


MRS 






3460000 


45 




1.010753 


MHMC 


.045 


26 


595000 


73 




1.010753 


MHMC 


.045 


26 


626000 


83 




1.010753 


MHMC 


.045 


20 


1044000 


189 


18 


1.010900 


MHMC 


.042 


28 


632000 


23 




1.010900 


MHMC 


.042 


28 


905700 


68 




1.010900 


MHMC 


.042 


28 


980000 


61 




1.010900 


MHMC 


.042 


28 


380000 


15 



measured probability distribution P{E) to differ- 
ent couplings j3 (see Fig p. To calculate the es- 
timates of our cumulants at each lattice size L 
we proceed in two steps: i) determine the error 
of each individual run performing a jackknife er- 
ror analysis by subdivision of the run into ten 
blocks; ii) calculate the final result by x^-fitting 
these individual run results to a constant. The 
L-dependence of the latter is quoted in Table ^ 
and|. 

4. FINITE SIZE SCALING ANALYSIS 

Under the assumption that the phase transition 
is discontinuous, the Borgs-Kotecky representa- 
tion of the partition function suggests, that both. 



the maxima of the specific heat 



and the 



pseudocritical /3- values can be expanded in terms 
of th in the inverse volume. 
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Figure 1. The measured cumulants of our data 
at L = 18 taken from the third run in table |^ 
reweighted to different couphngs /? in order to 
determine their extrema. 



and 

oo 

Pc{V)^l3T + Y.BkV-\ (8) 
fe=i 

The quantity G stands for the infinite volume gap 
in the plaquette energy and (3t denotes the infi- 
nite volume transition point of the system. We 
fit our data for C^^max and /3c (given in Table ||) 
to the parameterizations, Eq. ^ and |[ 

In order to expose systematic effects in the fit 
parameter G, we vary the fit range with various 
ranges Lmin < i < 18 as well as the truncation 
parameter kmax- The results in G are given in 
Table ^. It is possible to apply the parameteriza- 
tion of Eq. down to Lmm = 6 if one increases 



Table 2 

Locations of the specific heat maxima as a func- 
tion of lattice size. 



L 


/3c 


C\,,r.aj6V X 10-* 


6 


1.001794(64) 


9.7728(497) 


8 


1.007413(14) 


5.5428(105) 


10 


1.009383(16) 


3.8535(130) 


12 


1.010229(14) 


2.9813(118) 


14 


1.010626(11) 


2.5470(139) 


16 


1.010840(9) 


2.2620(161) 


18 


1.010943(8) 


2.1154(152) 



Table 3 

Locations of the Ui minima. 



L 




U2,^in X 10 


6 


1.001383(63) 


-2.5815(136) 


8 


1.007306(14) 


-1.3954(27) 


10 


1.009344(16) 


-0.9481(33) 


12 


1.010212(14) 


-0.7241(29) 


14 


1.010617(11) 


-0.6142(35) 


16 


1.010836(9) 


-0.5428(40) 


18 


1.010940(8) 


-0.5063(37) 


L 


f3u. 


UAjinin X 10 


6 


1.001174(63) 


-3.4489(182) 


8 


1.007241(14) 


-1.8627(36) 


10 


1.009319(16) 


-1.2651(44) 


12 


1.010200(14) 


-0.9661(39) 


14 


1.010610(11) 


-0.8194(46) 


16 


1.010832(9) 


-0.7241(53) 


18 


1.010938(8) 


-0.6753(50) 



kmax such that one degree of freedom remains. 
It is nice to see that results with Xdof < 1 
completely consistent. 

To illustrate the fit stability we additionally in- 
cluded the lowest expansion coefhcients Ai into 
Table |[ Note that here again Xdof < 1 results 
exhibit complete stability. 

The observed stability pattern supports the va- 
lidity of the ^-expansion. This corroborates 
the previous analysis R] on lattices up to L = 12. 
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Table 4 

Fit values of the infinite volume gap G according 
to Eq. 0. Fit results in bold face correspond to 
good quality fits (x^^f < !)■ 



Range 


^max 


Xdof 


G 


Ai 


14-18 


1 


1.03 


0.02731(17) 


2.62(12) 


12-18 


1 


6.65 


02779('27~) 


2.20(11) 




2 


0.50 


0.02698(18) 


3.22(20) 


10-18 


1 


25.9 


0.02837(37) 


1.88(9) 




2 


7.76 


0.02746(22) 


2.61(15) 




3 


0.72 


0.02686(26) 


3.47(34) 


8-18 


1 


202 


0.02965(73) 


1.40(8) 




2 


7.93 


0.02788(24) 


2.26(9) 




3 


2.21 


0.02737(22) 


2.72(16) 




4 


0.80 


0.02682(29) 


3.57(41) 


6-18 


1 


681 


0.03066(12) 


1.22(12) 




2 


113 


0.02913(60) 


1.64(10) 




3 


5.54 


0.02777(21) 


2.36(9) 




4 


2.17 


0.02735(22) 


2.75(17) 




5 


0.83 


0.02681(31) 


3.60(44) 



Table 5 

Infinite volume transition coupling (3t obtained 
by fitting Eq. ^. Fit results in bold face corre- 
spond to good quality fits. 



Range 


h 


Xdof 


Pt 


Bi 


14-18 


1 


0.47 


1.0111220(95) 


-18.81(56) 


12-18 


1 


0.36 


1.0111169(56) 


-18.46(24) 




2 


0.59 


1.0111238(166) 


-19.14(151) 


10-18 


1 


3.62 


1.0110978(124) 


-17.39(33) 




2 


0.83 


1.0111288(123) 


-19.65(44) 




3 


0.67 


1.0111207(196) 


-18.71(210) 


8-18 


1 


45.5 


1.0110437(370) 


-15.13(43) 




2 


0.44 


1.0111212(80) 


-19.08(20) 




3 


0.36 


1.0111299(82) 


-19.79(57) 


6-18 


1 


224 


1.0109952(785) 


-13.99(78) 




2 


17.8 


1.0110760(244) 


-16.72(42) 




3 


0.29 


1.0111257(45) 


-19.44(18) 




4 


0.36 


1.0111264(85) 


-19.81(62) 



As our final result for the infinite volume gap 
we obtain 



G = 0.02698(18)(17). 



(9) 



We proceed similarly in fitting the pseudocritical 
coupling /3t, Eq. ||, with results given in Table 
1^. From our analysis of Gy we obtain the critical 
Wilson coupling /3t = 1.011122(10)(8). 

As a further consistency check we analyze in 
addition the cumulants U2 and with param- 
eterizations analogous to Eq. ^. All three esti- 
mates (see Table ||) are in complete agreement 
within their statistical errors. We view this as a 
confirmation of the validity of Borgs-Kotecky FSS 
as applied to compact QED. As final we quote the 
average of the values given in Table @. 



Pt = 1.011128(11) 



(10) 



It remains to be seen whether our data exclude 
the possibility of a continuous singular behaviour 
with a vanishing infinite volume gap. Under this 
alternative assumption one expects the asymp- 
totic scaling laws 



Table 6 

Results for the transition coupling (3t obtained 
by analysis of three different cumulants. 



cumulant 


Pt 


Cv 


1.011122(10)(8) 


U2 


1.011129(14)(8) 




1.011132(6)(10) 



and 

I3,{L) = /3t + aL- 



(12) 



Here v denotes the critical exponent of the corre- 
lation length as hyperscaling is assumed. 

Our data clearly disfavour the validity of Eq. 
|ll|for any range of L- values. For Eq. |l^ however 
we find the fit parameters as quoted in Table ^ 

Note that fit stability cannot be reached un- 
der variation of the fit intervals Lmin ^ L < \S, 
neither in /3t nor in u. However it cannot be ex- 
cluded, that corrections to scaling could lead to 
stable fits also in this scenario. 



Gv,max{L) — CiL" 



(11) 
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Table 7 

Fit results obtained from second-order scaling in 
Eq. 



Lmin 


Xdof 


Pt 




6 


5.12 


1.011247(26) 


3.20(5) 


8 


2.17 


1.011204(22) 


3.34(7) 


10 


0.69 


1.011158(20) 


3.61(11) 


12 


0.64 


1.011127(32) 


3.90(29) 



5. DIRECT APPROACH TO LATENT 
HEAT 

Having determined Pt very accurately within 
Borgs-Kotecky FSS theory we are now in the posi- 
tion to make a 'direct' measurement of the latent 
heat based on our determination of (3t ■ Using 
the conventional canonical algorithm we generate 
O(IO^) configurations of each metastable phase 
simulating systems of size L = 6, 8, . . . , 24, 28, 32 
at Pt as given in Eq. |l^ and locate the positions 
of the energy peaks in the confined and Coulomb 
phase denoted by Ei and £'2- As shown in Fig. 
H we fit our data for each phase separately to a 
Gaussian of the form aexp{—b{E — Eif). Note 
that one cannot extract the relative phase weight 
from Fig. ||, because it depicts the probability dis- 
tributions Pi{E) of two independent runs j = 1, 2 
without any flip. Fig. || shows these peak posi- 
tions Ei{L) as a function of lattice size L. It is a 
prediction of the Borgs-Kotecky scheme that they 



L=32 , 







0.62 0.625 0.63 0.635 0.64 0.645 0.65 0.655 0.66 
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Figure 2. Separate fit of each measured phase 
peak of the L = 32 system at /3t to a Gaussian. 



0.665 




E 0.65 - 



0.645 - 

0.64 - 

0.635 - 

0.63 - 




Figure 3. Peak positions Ei as a function of lat- 
tice size L. 



deviate from their infinite volume value Ei{oo) by 
exponentially small terms only. Accordingly we 
fit both branches to 

E,{L) ^ E,{oo) + a^e-''^^ , (13) 

achieving high fit quality with Xdof i ~ O-'iO and 
Xdof 2 = 0.58. Fig. ^ illustrates the data and 
the fitted branches on an expanded energy scale. 
Due to the wider distributed peak of the confine- 
ment phase the errors of the lower branch Ei{L) 
are larger than in case of E2{L). With the fit 
results Ei{oo) we are in the position to present a 
high precision measurement of the infinite volume 
energy gap, 

G = £;2(oo) -£;i(oo) = 0.02667(20). (14) 

The estimate obtained from first-order scaling 
of the specific heat as given in Eq. ^ is in agree- 
ment with this result. This presents a further sup- 
port to the applicability of Borgs-Kotecky FSS to 
U(l) gauge theory. 

6. SUMMARY AND CONCLUSIONS 

All cumulants investigated in our high statis- 
tics analysis at L = 6, 8, 10, 12, 14, 16, 18 are in 
accord with Borgs Kotecky first-order FSS. Our 
data does not favour second-order FSS as indi- 
cated by large and lacking convergence of the 
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Figure 4. Data fitted to Eq. ^ on an expanded 
energy scale. 
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critical exponent v. Within the framework of 
Borgs-Kotecky we determined the infinite volume 
transition coupling with a relative error of the or- 
der of 10~^. A direct investigation of the latent 
heat at Pt for system sizes up to L = 32 yields a 
non-vanishing infinite volume energy gap G with 
an error of 1%. The systematic error of the energy 
gap due to the error of Pt is of the same order as 
its statistical error. The gap G = 0.02698(18) (17) 
extracted from first-order FSS of the specific heat 
is in agreement with the direct measurement of G. 

The material presented here is a progress re- 
port. In a forthcoming paper we shall elaborate 
further on systematic effects. 



